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Emerging devastating diseases, such as Huanglongbing (HLB) and citrus canker, have 
caused tremendous losses to the citrus industry worldwide. Genetic engineering is a 
powerful approach that could allow us to increase citrus resistance against these diseases. 
The key to the success of this approach relies on a thorough understanding of defense 
mechanisms of citrus. Studies of Arabidopsis and other plants have provided a framework 
for us to better understand defense mechanisms of citrus. Salicylic acid (SA) is a key 
signaling molecule involved in basal defense and resistance (R) gene-mediated defense 
against broad-spectrum pathogens. The Arabidopsis gene NDR1 (NON-RACE-SPECIFIC 
DISEASE RESISTANCE 1) is a positive regulator of SA accumulation and is specifically 
required for signaling mediated by a subset of R genes upon recognition of their cognate 
pathogen effectors. Our bioinformatic analysis identified an ortholog of NDR1 from citrus, 
CsNDRI. Overexpression of CsNDRI complemented susceptibility conferred by the 
Arabidopsis ndr1-1 mutant to Pseudomonas syringae strains and also led to enhanced 
resistance to an oomycete pathogen Hyaloperonospora arabidopsidis. Such heightened 
resistance is associated with increased SA production and expression of the defense 
marker gene PATHOGENESIS RELATED 1 (PR1). In addition, we found that expression 
of PR1 and accumulation of SA were induced to modest levels in citrus infected with 
Candidatus Liberibacter asiaticus, the bacterial pathogen associated with HLB disease. 
Thus, our data suggest that CsNDRI is a functional ortholog of Arabidopsis NDR1. Since 
Ca. L. asiaticus infection only activates modest levels of defense responses in citrus, we 
propose that genetically increasing SA/NDR1-mediated pathways could potentially lead to 
enhanced resistance against HLB, citrus canker, and other destructive diseases challenging 
global citrus production. 
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INTRODUCTION 

Huanglongbing (HLB; also called citrus greening disease), citrus 
canker, and other emerging diseases have imposed serious threats 
to the citrus industry worldwide (Bove, 2006; Gottwald, 2007). Cit- 
rus canker is caused by the gram negative bacterium Xanthomonas 
axonopodis pv. citri. Wind and rain facilitate the dispersal of the 
pathogen and spreading of the disease. More than 16 million cit- 
rus trees have been destroyed in Florida in an effort to restrict 
the disease (Gaskalla, 2006). HLB is even more devastating than 
citrus canker since it is highly contagious and lethal to afflicted 
plants (Bove, 2006; Brlansky and Rogers, 2007; Gottwald, 2007). 
The parasitic bacterium Candidatus Liberibacter that lives in the 
phloem tissue of a citrus tree is believed to be the associating 
agent of HLB. The disease is transmitted by a small insect vector 
of the family Psyllidae. The growth of psyllids cannot yet be reli- 
ably controlled by conventional insecticide application (Ichinose 
etal., 2010). Without successful measures to control the causal 
pathogen and its transmission vector, HLB is endemic to a variety 
of citrus species and related plants. Therefore, it is imperative to 



develop strategies to contain and eventually eradicate HLB and 
other diseases challenging the production of citrus worldwide. 

Successful manipulation of citrus disease resistance relies on 
a thorough understanding of defense mechanisms of the plant. 
Although not well understood yet in citrus, defense mechanisms 
are well studied in other plants, in particular the model plant 
Arabidopsis thaliana, and have been shown to be relatively con- 
served among plants (Nishimura and Dangl, 2010). Therefore, 
prior knowledge of defense mechanisms from other plants should 
help us to better understand citrus defense and ultimately develop 
effective strategies to combat devastating diseases challenging the 
citrus industry. 

Plant defense can be preformed or induced. The preformed 
defense includes some existing physical structures and chemical 
compounds produced by plants before infection that can restrict 
pathogen invasion. The induced defense can be activated upon 
pathogen invasion, involving sophisticated surveillance systems 
to recognize general elicitors from pathogens and subsequently 
to activate basal defense. Much stronger defense can be induced 
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when host resistance (R) proteins specifically recognize their 
cognate pathogen effectors; thus such defense is also termed as 
R protein-mediated defense. The largest class of R proteins is rep- 
resented by a family of proteins that have two conserved domains, 
nucleotide binding site (NBS) and leucine-rich repeat (LRR; Mar- 
tin etal, 2003; McDowell and Simon, 2006). This class of R 
proteins can be further divided into two groups according to 
the N-terminal sequence, coiled-coil (CC)-NBS-LRR and Toll- 
interleukin-1 receptor (TIR)-NBS-LRR. Some CC-NBS-LRR 
type proteins are found to signal through NON-RACE-SPECIFIC 
DISEASE RESISTANCE 1 (NDR1). For instance, an ndrl mutation 
compromises resistance conferred by the CC-NBS-LRR proteins 
RPS2, RPM1, or RPS5 to Pseudomonas syringae expressing the 
avirulence effectors avrRpt2, avrB and avrRpml, or avrPph3, 
respectively (Century etal., 1995; Aarts etal, 1998). On the other 
hand, some TIR-NBS-LRR type proteins functionally require 
ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1). For instance, 
an edsl mutant is immune-compromised to P. syringae avrRpsA, 
resistance to which is conferred by the TIR-NBS-LRR protein 
RPS4 but not by TIR-NBS-LRR type R proteins (Aarts et al., 1998; 
Falk etal., 1999). These observations suggest a general rule that 
these two subgroups of R proteins can activate distinct down- 
stream signaling pathways. However, exceptions to this rule also 
exist for some other NBS-LRR R proteins (McDowell et al, 2000; 
Bittner-Eddy and Beynon, 2001; Xiao et al, 2001). 

The small phenolic molecule salicylic acid (SA) plays a key 
role in signaling both basal and R protein-mediated defense 
( Hammond- Kosack and Jones, 1996; Tsuda etal., 2008) and is 
involved in resistance against diverse pathogens and in response to 
various stress conditions (Malamy etal, 1990; Rassmussen etal., 
1991; Sharma and Davis, 1997; Tsuda et al, 2008). While increased 
SA accumulation and/or signaling lead to enhanced disease resis- 
tance, disrupting these processes by gene mutations or transgene 
expression result in compromised defense against pathogens (Dur- 
rant and Dong, 2004; Lu, 2009). Genes involved in SA-mediated 
defense can affect SA biosynthesis, accumulation, and/or signaling 
(Lu, 2009). For instance, SA INDUCTION-DEFICIENT 21 EDS 16, 
encodes isochorismate synthase contributing to the majority of SA 
biosynthesis (Wildermuth etal, 2001). Both NDR1 and EDS1 are 
known to act upstream of SA to regulate SA accumulation (Falk 
etal., 1999; Shapiro and Zhang, 2001). Downstream of SA signal- 
ing, NONEXPRESSOR OF PR GENES 1 (NPR1) acts as a signal 
transducer that regulates systemic acquired resistance, a long- 
lasting defense against broad- spectrum pathogens at the whole 
plant level (Cao etal, 1997; Ryals etal, 1997; Shah etal., 1997; 
Dong, 2004). Overexpression of NDR1, EDS1, NPR1, or several 
other SA regulators confers enhanced disease resistance to a range 
of pathogens in Arabidopsis and/or in other plants (Chern etal., 
2001; Fitzgerald et al, 2004; Lin et al., 2004; Makandar et al, 2006; 
Malnoy etal, 2007; Pegadaraju etal, 2007; Sandhu etal, 2009; 
Gao et al., 2010). Therefore, manipulation of SA-mediated defense 
has the potential to introduce broad- spectrum disease resistance 
in plants. 

NDR1 encodes a glycosyl-phosphatidyl inositol-anchored 
plasma membrane protein that belongs to a large protein fam- 
ily (Dormann et al, 1995; Varet et al, 2002; Coppinger et al., 2004; 
Zheng et al, 2004). A recent study implicates the function of NDR1 



in mediating plasma membrane-cell wall adhesion (Knepper et al., 
201 1). NDR 1 -like genes widely exist in different plants (Lee et al., 
2006; Chong etal., 2008; Cacas etal, 2011). Besides NDR1, some 
Arabidopsis homologs of NDR1 were shown to be highly induced 
by pathogen infection and/or to confer enhanced disease resis- 
tance to P. syringae when overexpressed (Varet etal., 2002, 2003; 
Coppinger etal, 2004; Zheng etal, 2004). Thus NDR1 and some 
members in the family are critical components of plant defense. 

In this study, we report the isolation and characterization of 
a functional ortholog of NDR1 in citrus, named CsNDRl. We 
found that overexpression of CsNDRl complements the suscep- 
tibility of Arabidopsis ndrl-1 mutant to P. syringae avrRpt2 and 
further confers enhanced disease resistance to P. syringae avr- 
Rps4, which normally is not affected by the endogenous NDRL 
Overexpression of CsNDRl also led to increased resistance to 
the oomycete pathogen Hyaloperonospora arabidopsidis (Hpa) iso- 
late Noco2. CsNDR 1 -induced disease resistance is associated with 
increased SA accumulation and expression of the defense marker 
gene PATHOGENESIS RELATED 1 (PR1) in the transgenic Ara- 
bidopsis plants. In addition, we found that citrus infected with 
Candidatus Liberibacter asiaticus, a pathogen associated with the 
HLB disease, expressed modestly increased CsNDRl and SA levels, 
compared with mock- treated plants. We propose that genetic engi- 
neering to enhance SA/NDR1 signaling pathway in citrus could 
potentially enhance its resistance to HLB, citrus canker, and other 
emerging diseases. 

MATERIALS AND METHODS 
PLANT MATERIALS 

Arabidopsis plants were grown in growth chambers with a 12 h 
light/ 12 h dark cycle, light intensity at 200 |xmol m -2 s _1 , 60% 
humidity, and 22° C. The ndrl-1 mutant was previously described 
(Century etal., 1995, 1997). Citrus plants, "Valencia" (Citrus sinen- 
sis [L.] Osbeck), were grown on the greenhouse bench and kept 
at 24° C under natural light conditions. Plants were irrigated as 
needed and fertilized every 3 weeks using a water-soluble fertil- 
izer mix, 20N-10P-20K (Peters Professional, The Scotts Company, 
Marysville, OH, USA). 

BI0INF0RMATIC ANALYSIS 

Basic Local Alignment Search Tools (BLAST) was used to search 
protein sequence databases for Arabidopsis 1 and Citrus sinen- 
sis 2 , using appropriate query sequences. Sequence alignment 
and phylogenetic analysis were performed with the MEGA pro- 
gram (version 5.05). To construct the phylogenetic tree, the 
neighbor-joining method with 1000 bootstrap replications was 
used. 

PATHOGEN INFECTIONS 

Pseudomonas syringae strains used in this study were previously 
described (Lee etal., 2008; Wang etal., 2011b). Bacteria were cul- 
tured at 28°C with Kings B medium (10 g proteose peptone, 1.5 g 
K 2 HP0 4 , 3.2 ml 1 M MgS0 4 , and 5 g glycerol/1) containing the 
appropriate antibiotics. Freshly cultured bacteria at the optical 



1 http ://blast. ncbi. nlm. nih. gov/Blast . cgi 
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density of 0.5-0.8 were harvested, washed once, and resuspended 
in 10 mM MgS04 to make the infection solution at the desired con- 
centrations. The fifth to seventh leaves of 25-day-old Arabidopsis 
plants were infected by infiltration with a 1-ml needleless syringe. 
For bacterial growth assay, six leaves selected from over 10 plants of 
each genotype were collected 3 days post- infiltration, bored with 
a core borer (6 mm in diameter), and ground for bacterial growth 
measurement as described previously (Lu etal, 2003). For hyper- 
sensitive response (HR) test, one-half of a leaf at the fifth, sixth, 
or seventh position was infiltrated with P. syringae pv. maculicola 
ES4326 avrRpt2 (Pma avrRpt2; OD 600 = 0.1) and scored 16-24 h 
post-infiltration for leaf collapse. Leaves infiltrated with 10 mM 
MgS04 or the isogenic virulent strain Pma (OD600 = 0.1) were 
used as controls. At least 16 leaves from different plants of each 
genotype were scored for the HR symptoms. The rate of HR was 
expressed by the percentile of the number of leaves that developed 
HR symptoms out of the total number of inoculated leaves. 

Hyaloperonospora arabidopsidis isolate Noco2 was a kind gift 
from S. Xiao at University of Maryland College Park. Strain prop- 
agation and preparation were conducted as previously described 
(Song etal., 2004; McDowell etal, 2010). Hpa Noco2 spores 
(5 x 10 4 spores/ml in water) were sprayed on 7-day-old soil- 
grown seedlings. Sporangiophores on both sides of cotyledons 
were counted 7 days post-inoculation. At least 50 cotyledons from 
each genotype were counted to derive the average number of 
sporangiophores per genotype. 

For Ca. L. asiaticus infection, 15 -month- old "Valencia" plants 
were inoculated by grafting with two bark- or bud-pieces and 
two leaf pieces from infected greenhouse-grown "Valencia" plants, 
which were tested PCR-positive for Ca. L. asiaticus and demon- 
strated symptoms for HLB. Plants similarly inoculated but with 
disease-free tissue pieces obtained from healthy greenhouse-grown 
"Valencia" plants were used as controls. Plants were pruned imme- 
diately after graft-inoculation to promote new leaf growth and 
HLB disease development. The inoculated plants were random- 
ized periodically on the greenhouse bench to minimize the effect 
of environment on their defense responses to Ca. L. asiaticus. 

PCR-DETECTI0N OF Ca. L asiaticus IN CITRUS 

Ca. L. asiaticus-infected "Valencia" plants began to show typical 
HLB symptoms, yellowing and blotchy mottling around 1 1 weeks 
after inoculation (wai), which progressed more severely later. The 
symptomatic leaves were collected at 11 and 16 wai and extracted 
for DNA followed by PCR-detection of Ca. L. asiaticus (Albrecht 
and Bowman, 2012). Specifically, 100 mg leaf tissue was ground 
for DNA extraction, using the Plant DNeasy Mini Kit (Qiagen, 
Valencia, CA, USA) according to the manufacturer's instructions. 
For detection of Ca. L. asiaticus, real-time PCR assays were per- 
formed using primers HLBas and HLBr and the probe HLBp as 
described (Li et al, 2006). Amplifications were performed over 40 
cycles using an ABI 7500 real-time PCR system (Applied Biosys- 
tems, Foster City, CA, USA) and the QuantiTect Probe PCR Kit 
(Qiagen) according to the manufacturer's instructions. All reac- 
tions were carried out in a 20-|xl reaction volume using 5 |xl DNA. 
Once a branch was confirmed positive for Ca. L. asiaticus, symp- 
tomatic leaves on the branch were harvested for further RNA and 
SA analyses. 



RNA EXTRACTION AND ANALYSIS 

Twenty- five -day- old Arabidopsis plants were harvested for RNA 
extraction and northern blotting as described (Ng etal., 2011). 
Radioactive probes were made by PCR with an antisense primer 
specific for a gene fragment in the presence of [32P] dCTP. 
Primers used for making the CsNDRl probe were CsNDRl - 
Fl (ATGTCAGAAAACGCCGGTG) and CsNDRl -Rl (TTAAG- 
CAAAAATCAAGACAAAAAAATAC). Primers for the PR1 and 
18S rRNA probes were described previously (Ng et al., 201 1). 

Fully expanded leaves from Ca. L. asiaticus-infected "Valencia" 
plants showing HLB symptoms were collected 16 wai. One gram 
leaf tissue was ground in liquid nitrogen with a mortar and pes- 
tle and resuspended in 10 ml guanidinium isothiocyanate buffer 
(Chomczynski and Sacci, 1987). Total RNA was extracted as previ- 
ously described (Strommer etal, 1993) with slight modifications. 
Phenol/ chloroform/isoamyl alcohol (25:24:1) extraction was fol- 
lowed by two extractions with chloroform/isoamyl alcohol and 
precipitation of RNA with isopropanol at — 20°C overnight. RNA 
was pelleted by centrifugation at 10,000 g and 4°C for 1 h, resus- 
pended in 5 ml water and precipitated overnight at 0°C with an 
equal volume of 8 M LiCl. After centrifugation at 10,000 g and 4°C 
for 1 h, RNA was washed twice with 70% ethanol, air-dried, and 
dissolved in 500 |xl of water. RNA was further purified, using the 
RNeasy® MinElute Cleanup kit (Qiagen) according to the man- 
ufacturer's instructions. Total RNA was DNase-treated using the 
TURBO DNA-/ree-Kit™ (Ambion, Austin, TX, USA) according 
to the manufacturer's instructions. 

Quantitative reverse transcription real-time PCR (qRT-PCR) 
was performed using an ABI 7500 real-time PCR system 
(Applied Biosystems) and the QuantiTect SYBR Green RT- 
PCR Kit (Qiagen) according to the manufacturer's instruc- 
tions. Sixty nanograms of DNase-treated RNA were used in 
a total volume of 20 |xl. For detection of CsNDRl tran- 
scripts, forward primer S'-TGCTGCAGCTTCATCTTCAC-S' 
and reverse primer S'-TGTCGTGTTGTTTCGGTTGT-S 7 were 
used. For detection of 18S rRNA transcripts, forward primer 
S'-GCTTAGGCCAAGGAAGTTTG^ and reverse primer 5'- 
TCTATCCCCATCACGATGAA-3' were used. Melting curve anal- 
ysis was performed to ensure amplification of a single product 
and the absence of primer-dimers. For relative quantification of 
gene expression, the 2 _AACT method was applied as previously 
described (Livak and Schmittgen, 2001), using cycle threshold (Ct) 
values of 18S rRNA for normalization. 

cDNA AMPLIFICATION, DNA CONSTRUCTION, AND PLANT 
TRANSFORMATION 

To obtain the full-length CsNDRl cDNA sequence, we used 
the SMARTer™ RACE cDNA Amplification Kit (Clontech) to 
make a cDNA library from RNA extracted from Ca. L. asiaticus- 
infected 15-month-old "Valencia" plants. Nested primers, NDR1- 
5R-P1 (CACTTTCTGATCGGTCAGCGCAG) and NDR1-5R- 
P2 (CAATCACGGACGTGCCGATG), were used to amplify 
the 5 r end missing sequence while NDR1-3R-P1 (CATCG- 
GCACGTCCGTGATTG) and NDR1-3R-P2 (CTGCGCTGACC- 
GATCAGAAAGTG) were used to amplify the 3 r end missing 
sequence. The amplified fragments were cloned in the pJET 
cloning vector, using the Clone JET™ PCR Cloning Kit (Thermo 
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Scientific), and sequenced to obtain the full-length cDNA 
sequence. The full-length CsNDRl cDNA was further amplified 
from the library with NDR1-F1 (ATGTCAGAAAACGCCGGTG) 
and NDR1-R1 (TTAAGCAAAAATCAAGACAAAAAAATAC) and 
cloned into the pJET vector. At least 10 individual colonies were 
prepared for DNA and analyzed by sequencing. The sequence with 
fewer polymorphisms, compared with the reference sequence, and 
a correct open reading frame was used as the template for further 
cloning into the binary vector pBINplusARS under the control 
of the CAMV 35s promoter. The construct was confirmed by 
sequencing and transferred to Agrobacterium tumefaciens for ndrl - 
1 transformation, using the floral dipping method (Clough and 
Bent, 1998). T 0 seeds were selected for Ti plants on MS plates 
containing kanamycin, resistance to which was conferred by the 
binary vector. Ti transgenic plants were collected for seeds, which 
were further selected for homozygous T2 plants. 

ION LEAKAGE MEASUREMENT 

Leaves of 25-day-old Arabidopsis plants were infiltrated with a 
bacterial suspension Pma avrRptl (ODsoo = 0.1) with a blunt- 
end syringe, using 10 mM MgS04 treatment as a control. At 0, 
4, and 7 h post-inoculation, five leaf disks, cut with a 6-mm core 
borer, were collected, washed in de-ionized water, and placed in 
a 15-ml tube with 5 ml of de-ionized water. Triplicate tubes for 
each sample were gently shaken for 15 min followed by the mea- 
surement of solution conductivity, using an electrical conductivity 
meter (The London Company, Welwyn International Inc. Cleve- 
land, OH, USA). Each tube was measured three times to derive 
average conductivity. 

SA MEASUREMENT 

Free and total SA (glucosylated SA) were extracted from 25- 
day-old Arabidopsis plants (Ng etal., 2011; Wang etal, 2011a). 
The same protocol was used to extract SA from leaves of Ca. 
L. asiaticus-positive "Valencia" plants that demonstrated HLB 
symptoms. SA separation and detection were conducted with a 
high-performance liquid chromatography (HPLC) instrument as 
previously described (Ng et al., 201 1; Wang et al., 201 la). 

RESULTS 

IDENTIFICATION OF A CITRUS NDR1 0RTH0L0G 

Since NDR1 plays a critical role in Arabidopsis defense, we set out 
to identify the citrus ortholog and investigate its role in defense 
regulation. In Arabidopsis, NDR1 belongs to a large protein family 
with over 40 members, named NDRl/HINl-like (NHL) pro- 
teins (Dormanna etal., 2000). BLAST searching of the sequence 
database of Citrus sinensis 3 with the NDR1 protein sequence 
revealed a citrus protein (CsNDRl; orange l.lg0287 12m) with the 
highest similarity to NDR1 (the E- value is 2.4e -53 ) and three 
other top hits with E- values below 1.0 e -5 . We further used 
CsNDRl as the query to search the Arabidopsis protein database 
and retrieved sequences of NDR1 and 14 NHL proteins as the 
top hits. To determine the extent of similarity among these pro- 
teins, phylogenetic analysis was conducted, using the MEGA 
program (version 5.0; Tamura etal., 2011). Figure 1 shows that 



3 http ://www.phytozome.net 



CsNDRl is in the same cluster with NDR1 with 99% bootstrap 
support. Two other citrus NHL proteins (orange 1.1 g04 1808m 
and orange l.lg087 13m) are also in the same cluster with NDR1 
but with lower confidence levels in bootstrap support. Thus, 
bioinformatic analysis suggests that CsNDRl is an ortholog 
of NDRL 

ECTOPIC EXPRESSION OF CsNDRl COMPLEMENTS Arabidopsis ndrl-1 
MUTANT 

To test if CsNDRl shares conserved function with its Arabidopsis 
correspondence, we used a genetic complementation approach. 
The full-length CsNDRl cDNA was amplified via RT-PCR from 
a cDNA library made from Ca. L. asiaticus-infected "Valencia" 
plants and was cloned initially to the pJET vector and then to 
the binary vector pBINplusARS under the control of the CAMV 
35s promoter. The CsNDRl/ pBINplusARS construct was used to 
transform ndrl-1 via the standard floral dipping method (Clough 
and Bent, 1998). The presence of the transgene was confirmed 
by PCR with gene-specific primers. Initial infection of the Ti 
transgenic plants with a virulent strain R syringae pv. maculi- 
cola ES4326 (Pma) indicated that some of the transgenic plants 
were more resistant than ndrl-1 (data not shown). We further 
isolated homozygous plants for eight independently transformed 
lines (ndrl-1 + CsNDRl). Infection of these plants with Pma 
showed that all lines were more resistant than ndrl-1 and some 
were even more resistant than Col-0 (Figure 2A). Total RNA was 
isolated from these plants and northern blotting indicated that 
the level of disease resistance in some transgenic plants was corre- 
lated with the degree of transgene expression (Figure 2B). Thus, 
our results suggest that CsNDRl positively regulates Arabidopsis 
defense. 

The Arabidopsis RPS2 is a CC-NBS-LRR type R protein. When 
recognizing the avirulent strain Pma avrRpt2, RPS2 activates 
strong defense responses. Such defense activation requires the 
function of NDR1 and sometimes leads to HR, a rapid pro- 
gramed cell death in the infected region (Aarts etal., 1998). 
We found that all transgenic plants showed enhanced disease 
resistance to Pma avrRpt2 (OD 60 o = 0.0004), compared with 
ndrl-1 (Figure 3A). In addition, the ndrl-1 mutant showed 
compromised HR in response to a high dose of Pma avrRpt2 
(ODsoo = 0.1), as indicated by the lack of leaf collapse (Figure 3B, 
top panel; Century et al., 1995). We found that all transgenic plants 
showed partial to full rescue of HR-defect of ndrl-1 (Figure 3B, 
bottom panel) . 

Interestingly, line 15 that showed the highest level of CsNDRl 
expression, had a low frequency of leaf collapse in the HR assay, 
albeit still more than the ndrl-1 mutant. We also noticed that 
this line is smaller than other lines (Figure 4A). Thus we sus- 
pected that the small leaf size might obscure the HR scoring. 
To better quantify the HR cell death, we performed ion leakage 
measurement with this line and another line (#9) that showed 
medium CsNDRl expression (Figure 2B). When challenged with 
Pma avrRpt2 (OD 60 o = 0.1), ndrl-1 had the lowest level of 
ion leakage (Figure 4B; Zhang etal., 2004), consistent with its 
HR-deficit. The ion leakage level was highest in line 15 and 
medium in line 9, compared with Col-0. Together disease resis- 
tance and HR assays suggest that CsNDRl functions similarly as 
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FIGURE 1 | Neighbor-joining phylogenetic tree to show relationship 
among NDR1 homologs from Arabidopsis and citrus. The protein 
sequences were retrieved from BLAST search of Citrus sinensis protein 
database with NDR1 sequence as a query or from BLAST search of 
Arabidopsis protein database with CsNDRI sequence as a query, using an 
E-value cutoff e _4 .The MEGA program (version 5.05) was used to construct 



the tree. Protein sequences of the indicated genes were aligned with the 
ClustalW method and the tree was generated with the neighbor-joining 
method, using 1000 bootstrap replications. Numbers on the tree indicate 
bootstrap support (values <50% not shown). Branch lengths were drawn to 
scale; size bar represents number of amino acid substitutions per site. RIN4 
protein sequence was used as an outgroup to root the tree. 



Arabidopsis NDR1 in both basal and resistance protein-mediated 
defense. 

ECTOPIC EXPRESSION OF CsNDRI LEADS TO ACTIVATION OF 
SA-MEDIATED DEFENSE AND BROAD-SPECTRUM DISEASE 
RESISTANCE 

Salicylic acid is a key signaling molecule regulating defense path- 
ways including basal defense, R gene-mediated resistance, and 
systemic acquired resistance (Durrant and Dong, 2004; Lu, 2009). 
To test whether SA-mediated defense is activated in the transgenic 
plants, we quantified SA levels. We found that line 15 but not 
line 9 accumulated much higher levels of both free and total SA 
(glucosylated SA; Figure 4C). Consistent with its high SA levels, 
line 15 also showed higher expression of the SA marker gene PR1 



(Figure 4D). These results suggest that overexpression of CsNDRI 
to a certain level activates SA signaling. 

To further investigate how overexpressing CsNDRI affects 
disease resistance, we challenged line 9 and 15 with additional 
P. syringae strains. Ectopic expression of CsNDRI complemented 
susceptibility conferred by ndrl-1 to the virulent strain P. syringae 
pv. tomato DC3000 (Pto; Figure 5A, left). Line 15 is also more 
resistant to the isogenic avirulent strain Pto avrRps4 (Figure 5 A, 
right), which is recognized by RPS4 (a TIR-NBS-LRR type of 
R protein) independently of NDR1 (Aarts etal, 1998). Thus, 
these results indicate that ectopic expression of CsNDRI leads to 
activation of resistance to a pathogen that the endogenous gene 
otherwise does not have an effect on. In addition, we found that 
line 9 and 15 showed increased resistance to the virulent oomycete 
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FIGURE 2 | CsNDRI rescues enhanced disease susceptibility conferred 
by Arabidopsis ndr1-1 to a virulent P. syringae strain. (A) Bacterial 
growth assay. Twenty-five-day-old plants were infected with the virulent 
strain R syringae pv. maculicola ES4326 (Pma; OD 60 n = 0.0001). Bacterial 
growth was assessed 3 days after infection. Data represent the average of 
bacterial numbers in six samples ± standard error. Statistical analysis was 
performed with log transformed data, using one-way analysis of variance 
(ANOVA) Fisher's protected least significant difference (PLSD) tests 
(StatView 5.0.1). Different letters indicate significant difference among the 
samples (P < 0.05). (B) Expression analysis of transgenic plants carrying 
CsNDRI. Total RNA was extracted from 25-day-old plants and analyzed by 
northern blotting. A probe specific to CsNDRI was used to detect 
expression of the transgene. The 18S rRNA probe was used to indicate 
equal loading in the samples. These experiments were repeated two 
times with similar results. 



pathogen Hpa Noco2, compared with ndrl-1 (Figure 5B). Thus, 
overexpressing CsNDRI could lead to broad- spectrum disease 
resistance. 

SA ACCUMULATION AND EXPRESSION OF CsNDRI ARE MODESTLY 
INDUCED BY Ca. L asiaticus INFECTION 

To see how SA signaling is affected by HLB in citrus, we infected 
15-month-old "Valencia" plants with Ca. L. asiaticus, using mock- 
treated plants as a control. We began to observe at 1 1 wai HLB 
symptoms, chlorosis and/or blotchy mottling of leaves, which 
increased in severity by 16 wai (Figure 6A). We did qPCR with 
the symptomatic and control leaves, using primers specific to Ca. 
L. asiaticus. The average Ct values of the symptomatic leaves were 
21.8 at 1 1 wai and 20.7 at 16 wai. No Ca. L. asiaticus was detected 
in the control. Thus these symptomatic leaves were confirmed to 
be Ca. L. asiaticus positive. The symptomatic leaves were further 
collected for SA measurement and RNA analysis. Compared with 
the mock control, the symptomatic leaves from infected plants 
had about twofold more total SA levels (Figure 6B). Similarly, 
expression of CsNDRI was also induced about twofold more in 
the infected leaves (Figure 6C). These data suggest that SA and/or 
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FIGURE 3 | CsNDRI rescues enhance disease susceptibility and 
HR-defect conferred by ndr1-1 in response to Pma avrRpt2. (A) Bacterial 
growth assay. Twenty-five-day-old plants were infected with Pma avrRpt2 
(OD 60 o = 0.0004). Bacterial growth was assessed 3 days after infection. 
Data represent the average of bacterial numbers in six samples ± standard 
error. Statistical analysis was performed with log transformed data, using 
one-way ANOVA Fisher's PLSD tests (StatView 5.0.1). Different letters 
indicate significant difference among the samples (P < 0.05). (B) HR test. 
Half leaves of the fourth to seventh leaves of 25-day-old plants were 
infiltrated with Pma avrRpt2 (OD 600 = 0.1) and scored for the HR 
symptoms 16 h post-infection. At least 16 leaves were infected in each 
sample. Top panel: picture of infected leaves. The arrowheads indicate 
leaf collapse, a typical HR symptom, in Col-0 and a representative 
complemented plant (line 9) but not in ndrl-1. Bottom panel: the 
percentage of the infected leaves of each genotype that showed HR. 
These experiments were repeated two times with similar results. 
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FIGURE 4 | CsNDRI overexpression leads to increased ion leakage, 
SA accumulation, and PR1 expression. (A) Picture of 25-day-old plants. 
Other transgenic lines resemble Col-0 (not shown). (B) Ion leakage 
measurement. The fourth to seventh leaves of 25-day-old plants were 
infiltrated with Pma avrRpt2 (OD 60 n = 0.1) and collected at the indicated 
times post-infection for ion leakage measurement. The data represent 
average of triplicate samples, each of which contained five individual leaf 
disks of 6 mm in diameter. (C) SA quantitation. Twenty-five-day-old plants 
were harvested for SA extraction followed by HPLC analysis. Asterisks 
indicate significant difference between the samples and Col-0 (P < 0.05; 
n = 3). (D) Northern blotting. DNA fragments specific to PR1 and 18s 
rRNA, respectively, were used to probe blots containing RNA samples from 
Col-0, ndrl-1, and transgenic lines 9 and 15. These experiments were 
repeated two times with similar results. 



CsNDRI -mediated signaling are activated modestly upon Ca. L. 
asiaticus infection. 

DISCUSSION 

In this study, we presented bioinformatic and experimental 
evidence suggesting that the citrus gene CsNDRI is an Arabidop- 
sis NDR1 ortholog. Overexpression of CsNDRI in Arabidopsis 
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FIGURE 5 | CsNDRI confers broad-spectrum disease resistance. 

(A) Bacterial growth assay. Twenty-five-day-old plants were infected 
with Pto (left panel; OD 600 = 0.0002) or Pto avrRps4 (right panel; 
OD 60 o = 0.0004). Bacterial growth was assessed 3 days after infection. 
Data represent the average of bacterial numbers in six samples ± standard 
error. (B) Hpa sporangiophore measurement. Seven-day-old seedlings were 
spray-infected with Hpa Noco2 (5 x 10 4 spores/ml in water). Statistical 
analysis was performed with one-way ANOVA Fisher's PLSD tests 
(StatView 5.0.1). Different letters indicate significant difference among the 
samples (P < 0.05). These experiments were repeated twice with similar 
results. 



rescues ndrl-1 -conferred susceptibility to P. syringae infection 
and leads to broad- spectrum disease resistance to the oomycete 
pathogen Hpa Noco2. We also found that both SA accumula- 
tion and CsNDRI expression are induced to modest levels in 
citrus upon infection with Ca. L. asiaticus, the agent associated 
with HLB. Our data suggest a possibility that manipulation of 
SA/ CsNDR 1 -mediated defense may lead to enhanced resistance to 
HLB and other devastating diseases in citrus. 

SA plays a critical role in regulating plant resistance against 
various pathogens. Broad-spectrum disease resistance has been 
successfully introduced into several economically important plants 
via manipulation of the SA pathway. Much of the previous studies 
have been focused on NPR1, a key SA signal transducer. For 
instance, overexpression of Arabidopsis NPR1 and/or its homologs 
from other plants confers resistance against diverse bacterial and 
fungal pathogens in Arabidopsis, apple, citrus, soybean, tomato, 
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FIGURE 6 | Ca. L asiaticus infection induces modest increase of SA 
accumulation and CsNDRI expression in citrus. Fifteen-month-old 
"Valencia" plants were infected with Ca. L. asiaticus or mock treated. 
(A) Photo of leaves collected at 16 wai. Note the infected leaves showing 
chlorosis and blotchy mottling, typical HLB symptoms. (B) SA quantitation. 
(C) qRT-PCR analysis. CsNDRI expression was normalized to the 
expression of 18S rRNA and is presented as mean fold differences relative 
to control. Results are based on six biological replications. Statistical 
analysis was performed with one-way ANOVA Fisher's PLSD tests 
(StatView 5.0.1). Asterisks indicate significant different between the 
samples (P < 0.05). These experiments were repeated two times with 
similar results. 



rice, and/or wheat (Fitzgerald et al., 2004; Li et al., 2004; Makandar 
etal, 2006; Malnoy etal, 2007; Sandhu etal, 2009; Zhang etal, 
2010). These observations suggest that SA-mediated defense is 
conserved in monocots and dicots, and can be activated to against 
a range of pathogens (Campbell et al., 2002). Consistent with this 
notion, reports showed that exogenous application of SA agonists, 
such as benzothiadiazole and its commercial forms, acibenzolar- 
S-methyl (ASM, Actigard®, Syngenta Crop Protection, Inc.) and 
imidacloprid (Imid, Admire®, Bayer Crop Science), to citrus and 
other plants could induce defense marker gene expression and/or 
activate some protections of plants to a variety of viral, bacterial, 
and fungal pathogens (Campbell etal., 2002; Maxson- Stein etal., 
2002; Dekkers et al, 2004; Francis et al, 2009; Graham and Myers, 
2009). 

Arabidopsis NDR1 is a positive SA regulator, which is known 
to be specifically required for defense activated by some CC- 
NBS-LRR R proteins but not by some TIR-NBS-LRR type of R 
proteins (Aarts etal., 1998; Coppinger etal., 2004). Although not 
much is known about R-avr recognition in many non- Arabidopsis 
plants, the fact that NDR1 homologs widely exist in diverse plants 



(Lee etal, 2006; Chong etal, 2008; Cacas etal, 2011) suggests 
conserved defense signaling involving NDR1 homologs. Recently 
a coffee NDR1 homolog was shown to complement the ndrl-1 
mutant for its susceptibility to R syringae strains (Cacas etal, 
2011). Here we also showed a complementation of ndrl-1 by 
CsNDRI. The transgenic plants showed varying levels of CsNDRI 
expression, which is not uncommon for transgenes (Lu etal., 
2003). We noticed that there is a degree of correlation between 
the level of transgene expression and the level of disease resistance 
in some transgenic plants (Figure 2). One line highly expressing 
CsNDRI (line 15) constitutively activates SA-mediated defense, 
associated with enhanced disease resistance to Pto avrRps4 (that 
is not affected by the endogenous Arabidopsis NDR1) and to 
the oomycete isolate Hpa Noco2 (Figure 5). It is known that 
signals induced by different R genes upon recognition of their cog- 
nate effectors from pathogens can converge at downstream steps, 
involving SA-mediated defense (Martin etal., 2003; Chisholm 
et al, 2006). Our data suggest that hyper- activation of one branch 
of R-gene pathways, such as the CsNDRI branch, could poten- 
tially activate SA signaling, leading to broad- spectrum disease 
resistance. 

In Arabidopsis, both resistant and susceptible responses to 
pathogen infection are characterized by elevated SA accumu- 
lation and defense gene induction but with differences in the 
speed and amplitude of the responses (Zhou etal, 1998; Maleck 
etal, 2000; Tao etal, 2003; Song etal., 2004). Compared with 
a resistance response in Arabidopsis, induction of SA levels in 
citrus infected with Ca. L. asiaticus is quite small (Figure 6B). 
In addition, our gene expression data (Figure 6C) and microar- 
ray analyses (Albrecht and Bowman, 2008, 2012) indicate that 
the spectrum and intensity of defense genes induced by Ca. L. 
asiaticus are also quite limited. These observations suggest that 
when infected by Ca. L. asiaticus, citrus plants do not activate 
considerable host defense. This can be explained with at least 
two possible reasons: (1) a lack of recognition of effector pro- 
teins from Ca. L. asiaticus by citrus; and/or (2) a suppression of 
host defense by Ca. L. asiaticus. Thus, the interaction between 
citrus and Ca. L. asiaticus can be viewed as a compatible inter- 
action, leading to disease symptom development in the host. 
The relatively low level of host defense in response to Ca. L. 
asiaticus also suggests a possibility that HLB resistance can be 
achieved if we could manipulate the host to enhance its defense 
levels. 

Genetic engineering is a particularly attractive approach to 
introduce disease resistance traits into citrus because citrus 
has long juvenile growth - it typically takes 5-15 years for 
a citrus plant to flower. In addition, most commercial citrus 
cultivars produce polyembryonic seeds asexually, which com- 
plicates the process of introducing novel traits into citrus via 
traditional breeding (Koltunow etal., 1996). Recently the Ara- 
bidopsis NPR1 gene was shown to increase resistance to the 
canker disease when overexpressed in citrus (Zhang etal., 2010). 
Thus, CsNDRI, citrus NPR1, and other citrus homologs of 
SA regulatory genes are ideal candidates that can be geneti- 
cally manipulated to increase their expression in order to test if 
these genes confer resistance to HLB in citrus. Engineering such 
genes could yield citrus plants with enhanced disease resistance 
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